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Abstract 

The statistical mechanics of a non-interacting polymer chain in the limit of a large number of 
monomers is considered when the total angular momentum, L, is fixed. The radius of gyration for a 
ring polymer in this situation is derived exactly in closed form by functional integration techniques. 
Even when L = the radius of gyration differs from that of a random walk by a prefactor of order 
unity. The dependence on L is discussed qualitatively and the large L limit can be understood by 
physical arguments. 

PACS numbers: 82.35.Lr, 82.37.-j, 82.80.Ms 
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The statistical properties of polymers have been the subject of intensive research for 
many decades However these efforts have been almost entirely confined to polymers 
in liquids or solids, while in contrast, their properties in a vacuum have received little 
attention. Until recently there were not clear physical realizations of such situations, but 



now with recent developments in the mass spectrometry of proteins 



long hydrocarbon molecules in interstellar media [3j, such systems are now of experimental 



interest. The author (4] recently considered such systems theoretically and by means of 
computer simulation, and the purpose of this paper is to present an exact derivation for the 
radius of gyration of a polymer in a vacuum with conserved angular momentum. 

In reality angular momentum is weakly broken by interaction with thermal electromag- 
netic radiation [4J] but it is still important to understand the case of conservation laws 
properly in order to understand these more complicated effects. Also intra-chain interac- 
tions are not considered, that is this is the case of an "ideal" chain What we will find 
is at first sight rather surprising, that a polymer chain with conserved total energy E, total 
linear momentum p t ot = 0, and total momentum L, has a radius of gyration that depends 
strongly on L, so even when L = 0, the radius of gyration differs significantly from that of 
an ideal chain without this restriction. 

The statistical mechanics of a general classical system of N particles with constant total 
energy, momentum, and angular momentum with coordinates {r^} and momenta {p^} has 
been considered previously by Laliena [5j. They were considered to be interacting via a 
general potential $. In the case under consideration here, this microcanonical formulation 
can be expressed in a canonical ensemble in the limit of large N, after which it can be 
converted to a functional integral. This can then be evaluated to obtain the radius of 
gyration as a function of angular momentum L. 

Laliena ^] has shown that the conservation of linear momentum does not effect an- 
swers obtained in the microcanonical ensemble with conservation of angular momentum 
enforced {(J. So we will write down the volume of phase space with L and energy E kept 
constant. 

W{E,L,N) =C J 6(E-K-$)6<- 3 \L- y £T i x Pl )6^\r cm ) (jld 3 r i( i 3 p^ (!) 

C is a constant here that involves iV and h and is of no consequence for the purposes 
here. K is the kinetic energy J^p?/2m, with m the mass of each monomer, and here we are 
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and the search for 
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taking them all to be equal. The center of mass r cm also must be conserved and is set to 0. 
We use the Fourier representation of the 5 functions to write this as 

W(E,L,N) oc J c d\e XE J e - x ^5^(L -J> x Pl )5^(r cm ) f ft ^ft) 

As shown by Lax Q], for most purposes, as discussed below, the contour of integration can 
be deformed in the complex A plane using the method of steepest descents. The three con- 
ditions, (a) that a saddle point exists, (b) that an observable not be of order exp(const.N), 
and (c) that there be no singularity in the observable in the neighborhood of the saddle 
point. Conditions (b) and (c) are first obtained through the canonical ensemble and then 
tested to see if they are satisfied. Condition (a), that a saddle point exists, is satisfied 
because we can find a relationship between the energy and the temperature. In the case 
of an athermal system, say of rigid links, this would just be that E cx 1//?, where (3 is the 
value of A at the saddle point. Condition (b) is satisfied for the quantity of interest here, 
the average radius of gyration. Condition (c) is also satisfied because we will see that the 
average radius of gyration is smooth function of the temperature for T for finite T = 1/13. 
Thus we can drop the integration over A and replace A by the inverse temperature, /3, and 
consider the partition function Z instead instead of the phase space volume integral (which 
is simply related to the entropy.) 

Z(/3,L,N) oc Jd 3 kJ e ik - L e^^ + *)e-*-E,^xp i5 (3) (rcm) ^J[ d ^ p \ (3) 
Integrating over the pi's we obtain 

TV 

Z{/3, L, N) oc J d 3 k^ L J e^ ' V^ f[ d 3 n = J d^^CiP, k) (4) 

i=l 

where / is the moment of inertia tensor for the particles 
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Ia-y = m^{r u ir u i5 ai - rV 7 i) (5) 
i=i 

where a and 7 label the coordinates (1,2,3), and the Einstein summation convention has 
been used for v. In the last equality of Eq. @]we have introduced the function £(/3, k). Note 
that this cannot depend on the direction of k but only its magnitude, if $ only involves 
isotropic central potentials. Therefore we can take k to be along the z axis: k = kz, and 
write 

2 N 

COM) = / e" 1 *-^^-^^) Ud'n (6) 

J i=l 



Because of the radial dependence of £ on k we can also perform the k angular integrals in 
Eq. HI rewriting the k integration in spherical coordinates, obtaining 

c r°° 

Z((3,L,N) = - ksm(kL)(((3,k)dk (7) 
L Jo 

where c is a constant that plays no role in the subsequent analysis. 

The potential is taken to be that of an ideal Gaussian chain with step length I and a ring 
topology. 

^o = ^^E|r m -r i )| 2 + |r JV -r 1 | 2 j (8) 

By the central limit theorem, many models of polymer chains will all give the same results 
for most quantities of interest, if the overall radius of gyration is <C N. 

In order to calculate the radius of gyration, one can add an additional potential with a 
parameter e 

N 

/3$ = / 5$ + e /^|r i | 2 (9) 

i=i 

so that the average radius of gyration can be written as 

^ = <^I>I 2 > = - M — _ do) 

2=1 



e=0 



The integration in Eq. [6] is Gaussian and we can now take the usual limit to turn this into 
a functional integral 

C(&Jfc) = /e-r( z #^)^(^ 2 W)+- 2 +il r l^5(3)( rcm ) 5r(s ) (ii) 



The functional integration in the x, y, and z directions decouple and the x and y functional 
integrals are identical. Each one of these three integrals is of the form of the partition 
function of a one dimensional quantum harmonic oscillators at finite temperature except for 
the restriction on the center of mass. If we consider the Euclidean time action for a quantum 
harmonic oscillators of mass M at inverse temperature f3 Q 

S = M + „ o V) dt (12) 

then the partition function 

Z = f e- s 5x(t) oc -5 (13) 

J W 2sinh(^) V ; 
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with periodic boundary conditions on the paths x(0) = x(/3 ). This is true because of the 
general formula relating the partition function to the Euclidean path integral over times 
ranging from to the inverse temperature The partition function, Eq. [T3J, can also 
be derived in a less elegant but more direct manner by writing all paths in terms of a 
Fourier expansion, which then decouples the integrals, and forms an infinite product over 
all modes. This latter approach is useful in the present application because we have the 
additional restriction on the path integral that the zero mode should not be integrated over 
as a consequence of the restriction on the center of mass. Using the Fourier decomposition 
approach, we can easily incorporate this restriction, by not including the zero mode in the 
product. This amounts to multiplying Eq. d3]by uj 

Rescaling variables so that the angular momentum L' = Ly/\2/{Nl\/wT) and using Eq. 
[TU1 gives 

Rg 2 1 /o°°(fc(-6 + k 2 + 6fccoth(fc))csch(fc) 2 sm(kL'))dk 
NP~ 36 J™ k 3 csch{k) 2 sm{kL')dk ^ ' 

Both the numerator and denominator can be computed in closed form using contour 

integration. This gives the final result 

R g 2 _ 2L' (3 + 7T 2 ) + L' (tt 2 - 6) cosh[LV] + 3 (L' 2 - 1) TrsinhjLV] 

iVT 2 ~ 367r(2L'7r + L'vrcoshfL'vr] - 3sinh[L'7r]) ^ ' 

A plot of this equation is displayed in Fig. [TJ 

For V = this reduces to R 2 /(Nl 2 ) = (1 + 15/7r 2 )/36 rs 0.07. For a ring without angular 
momentum conservation, R 2 /(Nl 2 ) = (1/12) w 0.083 which means that the restriction to 
L = causes the rings to be smaller relative to the case where the angular momentum can 
take on any value. 

In the opposite limit of large L', R 2 /(Nl 2 ) — > L'/12ir. Note that in terms of L, R 2 — > 



LI j (it \jYXTm) independent of chain length N. To understand this behavior, we consider 
high L configurations, where we expect that a typical configuration of the ring will be close 
to a circle rotating rapidly. The approximate free energy contains a kinetic energy and an 
elastic term 

J? h 

F = — + V (16) 

Where C = 2nR is the circumference, and k is the entropic elastic spring coefficient k = 
3T/(Nl 2 ). The moment of inertia is approximately I = mNR 2 . Minimizing with respect to 
R 2 this gives the above result. It is not surprising that this result is exact because in the 
large L limit, we expect that this circular configuration will become dominant. 
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FIG. 1: The radius of gyration versus rescaled angular momentum for an ideal ring chain. 



The temperature in this model at a given energy is determined in the usual way, by 
requiring that the average energy in the canonical ensemble is equal to the microcanonical 
energy. However in this case, the radius of gyration is a very sensitive function of L. To 
see this, note that changes occur on a scale L' ~ 1, or L ~ NlVmT. Because / is typically 
~ (Nm)R 2 g ~ (iVm)(iV/ 2 ) ~ N 2 l 2 m then L ~ NlVrtT. The order of L 2 /2I is therefore ~ T. 
This means that a change of order one degree of freedom changes R 2 /(Nl 2 ) by a number of 
order unity, which has a negligible effect on the temperature but a large effect on the radius 
of gyration. So for fixed V as iV — > oo, we see that in the canonical ensemble, the effect of 
the angular momentum constraint on the energy is a fraction of order 1/N. This means that 
when the limit N — * oo is taken with V fixed, the relation between the temperature and 
energy can be obtained as it would for a polymer without angular momentum conservation, 
and thus will not have any dependence on V . This is also seen more rigorously by computing 
the exact dependence of the partition function on L', which is done below. 

It is useful to calculate the probability density for finding the polymer with a particular 
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value of total (rescaled) angular momentum V . This should be important in the case where 
there is a dilute gas of such polymers. It is also important for a single chain for long times, 
since the angular momentum is changed by the weak coupling to electromagnetic black 
body radiation [4]. In this case the probability density function P(L') is proportional to the 
partition function Z(f3,L',N). The normalization requirement is that 

f° P(L')AiTL' 2 dL' = 1 (17) 
Jo 

The normalization is straightforward to calculate using Eq. [7] and integrating over V 
first. The V integration requires evaluating 

r u sm(kL')dL' = - r dc °<kL') dL , = _ n5 ,( k) (18) 

Jo Jo dk 

and using this, the integral over k is now easily accomplished 

/;z (M 4.L^ OT ?<if% o (19) 

Z was evaluated previously in the process of calculating the radius of gyration and is 
proportional to 

poo 

/ fc 3 csch(fc) 2 sin(H7)dA;/L' (20) 
Jo 

Evaluating this integral and including the correct normalization using Eq. [19] yields 



. 7r 3 csch(^) 4 (2L'7r + L'it cosh(L / 7r) - 3 sinh(L'7r)) 
P{L ] ~ 16L^ (21) 



ln(P(L')) versus V is plotted in Fig. [2J Because of the non-constant value of the moment 
of inertia for the chain, this distribution is decidedly non-gaussian. In the large V limit, the 
slope of this curve approaches a constant with a slope of — n. This is in agreement with the 
minimization argument for large V given under Eq. [161 

It is interesting to compare the extreme sensitivity of this system to restrictions in angular 
momentum, with what would be expected in other kinds of systems. The system consid- 
ered here is essentially one dimensional in that interactions are only from nearest neighbor 
monomers. In, for example, a membrane or a three dimensional gel, the system is of higher 
dimension. In such two or three dimensional systems, a perturbation that changes the free 
energy by O(ksT) is expected to only effect averages by microscopic amounts. However for 
polymers, it has a much larger effect. For example, a force pulling the ends of a polymer 
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FIG. 2: The probability density function for finding a chain in thermal equilibrium with rescaled 
angular momentum L' for an ideal ring chain. 



costing k B T of energy will increase its radius of gyration by a multiplicative constant of 
order unity. Because a polymer chain with the restriction L = 0, reduces the number of 
degrees of freedom by at least one, it should affect the free energy by 0{ksT). So by the 
above argument, this is expected to make non-trivial changes to its statistics, unlike higher 
dimensional systems. So even in the case L = 0, the non- interacting polymer is no longer 
an uncorrelated random walks. 

The author wishes to thank Peter Young for useful discussions. 
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